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The spatially resolved broad-band spectroscopy with Very Long Baseline Interferometry (VLBI) is one of the 
few methods that can probe the physical conditions inside blazar jets. We report on measurements of the 
magnetic field strength in parsec-scale radio structures of selected bright Fermi blazars, based on fitting the 
synchrotron spectrum to VLBA images made at seven frequencies in a 4.6-43.2 GHz range. Upper limits of 
B± < 10~^— 10^ G (observer's frame) could be placed on the magnetic field strength in 13 sources. Hard radio 
spectra (—0.5 < a < -1-0.1, Si, ~ v") observed above the synchrotron peak may either be an indication of a hard 
energy spectrum of the emitting electron population or result from a significant inhomogeneity of the emitting 
region. 



1. Introduction 

Blazars are active galactic nuclei characterized 
by highly variable non-thermal continuum emission 
across the electromagnetic spectrum. The Large Area 
Telescope (LAT, fl|) on board the Fermi 7-ray obser- 
vatory, launched on 2008 June 11, provides a wealth 
of new information about high-energy radiation of 
blazars. It is believed that this high-energy emission 
is intimately connected with the emission at lower 
energies (radio - optical) through the process of in- 
verse Compton scattering. Very Long Baseline Inter- 
ferometry (VLBI), with its spectacular angular reso- 
lution, provides the most detailed view of inner re- 
gions in blazars ^J. In most cases, VLBI observations 
reveal one-sided relativistically beamed parsec-scale 
jets originating in a bright, compact region called the 
"core" . 

The spatially resolved broad-band spectroscopy is 
one of the few methods that can probe the physical 
conditions inside blazar jets. Despite the great po- 
tential of this method, its applications are still rela- 
tively rare (see |3i] and references therein), owing to 
difficulties in implementing it and interpreting results 
obtained 'i'l . In addition, multifrequency VLBI obser- 
vations require large amounts of observing time and 
most VLBI arrays lack technical capability to perform 
simultaneous multifrequency observations. 

We utilize the Very Long Baseline Array (VLBA, 
0) which has a unique capability to conduct ob- 
servations quasi-simultaneously at many frequencies, 
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needed to constrain physical conditions in the regions 
where blazar radio emission originates. We recon- 
struct radio spectra of parsec-scale features in jets of 
selected 7-ray-bright blazars and compare them with 
the simple homogeneous synchrotron source model, 
which allows us to derive information about the mag- 
netic field strength and particle energy distribution. 
During the first year of the Fermi scientific operations, 
we used the VLBA to observe 20 blazars, that were 
identified prior to the launch of Fermi as expected 
bright 7-ray emitters. In this paper, we present results 
for 17 of them. The data analysis for three remaining 
sources is still ongoing. 



2. Observations and data reduction 

The observations were conducted simultaneously at 
seven frequencies in a 4.6-43.2 GHz range. After the 
initial calibration in the package [Q], the sources 

were self-calibrated and imaged independently at each 
frequency using the Difmap software 0. Whenever 
feasible, the frequency-dependent core position (the 
"core shift" effect, see for a detailed discussion [1,[1]) 
was taken into account. Image alignment is a ma- 
jor problem for spectral VLBI analysis because the 
self-calibration technique used to obtain high-quality 
images causes loss of information about the absolute 
position of a source. Positions of optically thin jet 
features or the core itself (if the "core shift" effect is 
negligible for a given source) are then used to align 
VLBI images obtained at different frequencies. 

In order to reconstruct the spectrum in a given 
pixel of an image ("pixel-based approach"), the 
sum of CLEAN components around this pixel was 
used. The CLEAN components were summed in- 
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side a radius R = 9 (lowest frequency) /3, where 
(lowest frequency) is the half-power beam width at 
the lowest observing frequency (4.6 GHz). This ap- 
proach has a major advantage over the direct use 
of flux densities from reconstructed images because 
it avoids artifacts introduced by convolution with a 
Gaussian beam while preserving the angular resolu- 
tion achieved. 



3. Results 

The radio spectra in the parsec-scale core regions 
(marked as shaded areas on Fig. [1]) of the observed 
blazars have diverse properties. Most of them are 
nearly flat, with signs of curvature at lower fre- 
quencies, which we interpret as synchrotron self- 
absorption. In one source (NRAO 530, Fig. [Tr)) . 
the spectral turnover is particularly prominent and 
well sampled by our observations. Cores in three 
sources, 0716+714 {Fig.^, OJ 287 (Fig. [ig, and 
3C 454.3 (Fig. I2dp . show inverted spectra which are 
slightly curved. The core spectrum of only one source 
(1510—089, Fig.[2c]) can be adequately described by a 
simple power law in the whole 4.6-43.2 GHz frequency 
range. In the following, we focus on sources that show 
a well defined spectral turnover, and thus can be used 
to determine the physical conditions in the emitting 
plasma. 

We fit a homogeneous synchrotron source model 
[l0| (see the Appendix for more details) to the ob- 
served spectra in order to obtain the magnetic field 
strength and particle energy distribution. The sources 
for which the data reduction is complete are presented 
in Table m For 13 sources, the model provided an ad- 
equate fit to the observed spectra (Fig. [T]) while four 
sources exhibit flat or inverted spectra (Fig. [2]) that 
cannot be described by the homogeneous synchrotron 
model. 



4. Discussion 

4.1. Model reliability 

The physical parameters presented in Table |T] cor- 
respond to typical values for emitting electrons inside 
a large area of a few milliarcseconds in size (shaded 
areas on Fig.[T]). In most cases, only an upper limit on 
the magnetic field strength could be placed, because 
most regions where the spectral turnover is detected 
remain unresolved. The results are consistent with es- 
timates obtained by a method based on model-fitting 
bright VLBI components with 2D Gaussian compo- 
nents (see Fig. |3] for results on BL Lacertae). The 
model-fitting based method is described in detail in 
[1] where it was applied to 3C 273 and the average 
magnetic field strength B± ~ 10~^ G was obtained 



for the 8 GHz core region which is comparable in size 
to the region probed in our research. This value is in a 
good agreement with our upper limit at B± < 0.016 G 
given in Table HI 



4.2. Electron energy spectrum 

The spectra above the synchrotron turnover are 
nearly flat or slightly inverted in core regions of all 
sources investigated in this work. This may result 
from blending of a few emitting regions with different 
peak frequencies (an example of such blending may be 
found in the core region of BL Lacertae, see the top 
right panel in Fig. [31 see also [Sj] for 3C 273). Alter- 
natively, a nearly flat spectrum may imply a hard en- 
ergy spectrum of the emitting electrons (as may be the 
case for the component Bl in the jet of BL Lacertae 
whose spectrum is presented in the lower right panel 
of Fig. 131 see also Table llll) . The hard electron spec- 
trum is difficult to explain by the conventional first- 
order Fermi acceleration mechanism. However, the 
hard spectrum can be produced by the second-order 
Fermi mechanism ( "stochastic acceleration" ) 1 1] . 



5. Summary and prospects 



We have placed upper limits of B± < 10^^-10^ G 
(observer's frame) on the magnetic field strength in 
13 out of 17 Fermi blazars (Table |T]). Spectra of four 
blazars could not be described by the simple homoge- 
neous synchrotron source model, and no constraints 
on the magnetic field and particle energy spectrum 
could be obtained for these sources. Hard spectra 
(—0.5 < a < +0.1, S„ ~ z^") observed in the blazar 
cores above the synchrotron peak may either indicate 
a hard energy spectrum of the relativistic electron 
population in the jet or result from significant inhomo- 
geneity (of the magnetic field and plasma parameters 
or just optical depth) across the emitting region. 

Since 7-ray emission in blazars is suggested to orig- 
inate from regions spatially close to the VLBI core 
fl^l, the estimates of the magnetic field strength 
and electron energy distribution presented in Ta- 
ble U could be used to constrain broad-band Spectral 
Energy Distribution (SED) models. Simultaneously 
with the VLBA observations described here. Swift X- 
ray/ultraviolet/optical target of opportunity observa- 
tions were performed. When combined with the Fermi 
7-ray data this will provide a unique data set for the 
SED modeling. Analysis of this data set is presently 
underway. 
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Table I: Physical parameters estimation 



lAU Name (B1950.0) Alternative Name 


Epoch 




B^'" [G] 


Comments 




0235+164 


AO 0235+16 


2008-09-02 


0.8 


< 0.89 






0528+134 




2008-10-02 


1.4 


< 0.12 






0716+714 




2009-02-05 






inverted spectrum a 


= 0.4, model N/A"^ 


0827+243 


OJ 248 


2009-04-09 


1.6 


< 3.4 






0851+202 


OJ 287 


2009-02-02 






inverted spectrum a 


= 0.7, model N/A'' 


1219+285 


W Com 


2009-05-14 


0.8 


< 108 


preliminary analysis 




1226+023 


3C 273 


2009-02-05 


2.0 


< 0.016 






1253-055 


3C 279 


2009-02-02 


1.4 


< 0.89 






1510-089 




2009-04-09 






inverted spectrum a 


= 0.2, model N/A'' 


1633+383 


4C 38.41 


2009-06-20 


1.0 


< 0.24 


preliminary analysis 




1652+398 


Mrk 501 


2009-05-14 


1.6 


< 3.5 


preliminary analysis 




1730-130 


NRAO 530 


2009-06-20 


1.6 


< 3.6 


preliminary analysis 




1959+650 




2008-10-23 


1.4 


< 17 






2155-304 




2008-09-05 


1.0 


< 26 






2200+420 


BL Lac 


2008-09-02 


1.0 


< 0.44 






2251+158 


3C 454.3 


2008-10-02 






inverted spectrum a 


= 0.8, model N/A'' 


2344+514 




2008-10-23 


1.2 


< 57 


43 GHz data not included in the fit 



p is the power law index in the electron energy distribution N{E) = NoE~^. 
Note, that for the optically thin part of the synchrotron spectrum p = 1 — 2a, 
where a is defined as ^ i^" ■ 

^ The estimates correspond to the region of parsec-scale radio core. 
The values are in the observer's frame. 

The homogeneous synchrotron source model is not applicable for this spectrum. 



Acknowledgments 

K. Sokolovsky is supported by the International 

Max-Planck Research School (IMPRS) for Astron- 
omy and Astrophysics at the universities of Bonn and 
Cologne. T. Savolainen is a research fellow of the 
Alexander von Humboldt foundation. T. Savolainen 
also acknowledges support by the Academy of Fin- 
land grant 120516. Y. Kovalev was supported by the 
return fellowship of Alexander von Humboldt foun- 
dation and the Russian Foundation for Basic Re- 
search (RFBR) grant 08-02-00545. This research has 
been partially funded by the Bundesministerium fiir 
Wirtschaft und Technologic under Deutsches Zentrum 
fiir Luft- und Raumfahrt grant number 5GORG808. 
This work is based on data obtained from the National 
Radio Astronomy Observatory's Very Long Baseline 
Array (VLBA), project BK150. The National Radio 
Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. This research 
has made use of NASA's Astrophysics Data Sys- 
tem; the NASA/IPAC Extragalactic Database (NED) 
which is operated by the JPL, Caltech, under contract 
with the NASA, the SIMBAD database, operated at 
CDS, Strasbourg, France. The authors are grateful to 
E. Ros for reviewing this manuscript. 



References 

1 Atwood W. B., et al., 2009, ApJ, 697, 1071 

2 Zensus J. A., 1997, ARA&A, 35, 607 

3 Savolainen T., Wiik K., Valtaoja E., Tornikoski M. 
2008, ASPC, 386, 451 

4 Lobanov, A. P. 1998, A&AS,132, 261 

5 Napier P. J., 1994, lAUS, 158, 117 

6 Fomalont E., 1981, NRAON, 3, 3 

7 Shepherd M.C., Pearson T.J. and Taylor G.B. 
1994, BAAS, 26, 987 

8 Kovalev, Y. Y., Lobanov, A. P., Pushkarev, A. B., 
& Zensus, J. A. 2008, A&A, 483, 759 

9 Lobanov, A. P. 1998, A&A, 330, 79 

10 Pacholczyk A.G. 1970, Radio astrophysics, 
San Francisco: Freeman 

11 Virtanen, J. J. P., & Vainio, R. 2005, ApJ, 621, 
313 

12 Kovalev, Y. Y., et al. 2009, ApJL, 696, L17 

13 Marscher A.P. 1983, ApJ, 264, 296 



eConf C091122 



4 



2009 Fermi Symposium, Washington, D.C., Nov. 2-5 



Appendix: synchrotron source model 
used in pixel-based approach 

Following [lo| . we consider a uniform cloud of rel- 
ativistic electrons in the external uniform magnetic 
field with component perpendicular to the line of 
sight. We assume that the electron energy distribu- 
tion is described by a power law: 

N{E) = NoE-P, 



J(z,p) = z5/2(i _ exp[-z-(P+4)/2]), 



Cl 



3e 



C3 



An rUeC^ 



6.27 X 10l^ 



1.87 X 10-2^ 



where Nq and p are constants. For this case, we can 
write the radiation transfer equation 



ds 



C5{P) = tCsF 



3p- 1\ /3p + 7\ fp + 7/3 



12 



12 



where Ii, is the specific intensity at the frequency v, 
ds is the differential element of path length, ji, is the 
emission coefficient (emissivity) , k^, is the absorption 
(extinction) coefficient. If and Ki, are constant 
along the line of sight, the solution of the radiation 
transfer equation is 



— f Ki,ds 

i^e f ds 



(1 



where L is the size of the electron cloud along the line 
of sight, and k^L = is the optical depth. 

Following [l3|, this solution can be written in the 
form: 



S{vi)J{ — ,p), 
^1 



(1) 



where 



V, = 2c^{NoLc,{p)f/^P+^^B^l+^^'^P+'\ (2) 



C5(p) ^_l/2 Vl 5/2 
Cfi(p) ^ 2ci 



(3) 



10 

y 



3p + 2\ /3p+ 10 



12 



12 



where e is the electron charge, me is the electron mass, 
and c is the speed of light in vacuum. 

The power law index p of the electron energy dis- 
tribution may be found from the spectral index a 
{Su ~ J/") in the optically thin part of the spectrum: 
p = 1 — 2a. The optical depth t„i corresponding to 
the observed spectral peak frequency Vm may be found 
from equation exp[T„i] = 1 -f ^^t^. 

The frequency vi corresponding to r^^ = 1 may be 
found from Vm, P and Tm'- — Tm^^^^^^ Vm- 

From equations (HJ, ([2]) and ([3]), we can derive the 
perpendicular component of the magnetic field: 



Si) ^<-'iS)' 



(4) 



Note, that if the size of an emitting region is smaller 
then the angular resolution, only a lower limit on the 
intensity and, correspondingly, an upper limit on the 
magnetic field B± can be derived. 
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Figure 1: VLBI images at 8.4 GHz and radio continuum spectra of sources for which the synchrotron self- 
absorption turnover frequency was successfuUy measured. The shaded area on the images marks the region 
where the synchrotron turnover was detected. The color of the shaded area represents the upper limit of the 
magnetic field strength {B±). The values of B± and p presented in Table U correspond to this area. On the 
spectrum plots, the blue curves represent the homogeneous synchrotron source model, while the green lines 
represent a simple power law model. Best fit parameters of both models are shown on each plot. The f sym- 
bol marks sources for which the analysis results are preliminary. 
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Figure 1: continued 
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Frequency (GHz) 
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Figure 2: VLB A spectra of sources whose spectra cannot be adequately described by the simple homoge- 
neous model. The spectra correspond to parsec-scale core regions. Green lines represent a power law fit. It is 
possible to fit a homogeneous synchrotron source model (blue curve) to the observed spectrum of 0716+714, 
however, a simple power law provides better (lower x^) fit- 
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Figure 3: Results of the model- fitting based analysis of BL Lacertae (2200+420). The contours represent 
the CLEAN image at 43 GHz, the beam size and the linear scale are indicated at the lower left panel. The 
source structure was approximated with Gaussian components (see Table |TT| . The fitting was conducted in 
the the uv-plane. The two inner components (CI its position is marked by the red arrow, C2 marked by 
the blue arrow) are unresolved while the third component (Bl marked by two black arrows) is resolved at 
43 GHz. Fig. [Tx] presents the spectrum of this source obtained in the course of pixel-based analysis. This 
spectrum corresponds to the combined spectrum (black) of components CI (red) and C2 (blue) presented 
in the upper right panel of the current figure. The peak frequency for Bl is not well constrained, however, 
following [l^ we can put an upper limit on the magnetic field strength in this component B < 0.06 G (ob- 
server's frame). 



Table II: BL Lacertae model component parameters 



Component 


Distance 






a 




[mas] 


[GHz] 


[Jy] 




CI (core) 





> 43 


> 1.02 




C2 


0.26 


6.4 


1.43 


-0.58 


Bl 


1.47 


< 4.1 


> 0.51 


-0.17 



z-'m — synchrotron self-absorption peak frequency. 

Sm — synchrotron self-absorption peak flux density. 

a — spectral index in the optically thin part of the spectrum. 
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